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X am repeating the proposal abstract here to provide the context for 
the final report: 

I propose to continue providing observers with basic data for interpreting 
spectra from stars, novas, supernovas, clusters, and galaxies. These data 
will include allowed and forbidden line lists, both laboratory and computed, 
for the first five to ten ions of all atoms and for all relevant diatomic 
molecules. I will eventually expand to all ions of the first thirty elements 
to treat far UV and X-ray spectra, and for envelope opacities. I also include 
triatomic molecules provided by other researchers. X have made CDs with 
Partridge and Schwenke's water data for work on M stars. The line data also 
serve as input to my model atmosphere and synthesis programs that generate 
energy distributions, photometry, limb darkening, and spectra that can be used 
for planning observations and for fitting observed spectra. The spectrum 
synthesis programs produce detailed plots with the lines identified. Grids 
of stellar spectra can be used for radial velocity-, rotation-, or abundance 
templates and for population synthesis. X am fitting spectra of bright stars 
to test the data and to produce atlases to guide observers. For each star 
the whole spectrum is computed from the UV to the far IR. The line data, 
opacities, models, spectra, and programs are freely distributed on CDs and on 
my Web site and represent a unique resource for many NASA programs. 

X am now in full production of new line lists for atoms [, except that 
I have not been getting paid since Hay 15, 2004, and except that X am going 
to run out of money for supplies and equipment] . X am computing all ions of 
all elements from H to Zn and the first 5 ions of all the heavier elements, 
about 800 ions. For each ion X treat as many as 61 even and 61 odd 
configurations, computing all energy levels and eigenvectors. The 
Hamiltonian is determined from a scaled-Hartree-Fock starting guess by least 
squares fitting the observed energy levels. The average energy of each 
configuration is used in computing scaled-Thomas-Fermi-Dirac wavefunctions 
for each configuration which in turn are used to compute allowed and forbidden 
transition integrals. These are multiplied into the LS allowed and forbidden 
transition arrays. The transition arrays are transformed to the observed 
coupling to yield the allowed and forbidden line lists. Results are put on 
the web as they are finished. [Provided X get funding,] there will be more 
than 500 million lines. I will then compare ion by ion, to all the laboratory 
and computed data in the literature and make up a working line list for 
spectrum synthesis and opacity calculations with the best available data. 

As the laboratory spectrum analyses are improved, X will redo the calculations 
with the new energy levels. 

My original plan when X started the new calculations was to run through 
all the atoms using my old Cray programs from the 1980 's that were limited to 
1100 x 1100 arrays in the Hamiltonian for each J. Then X would go back and 
rerun the more complicated cases with 3000 x 3000 arrays so that X could 
include many more configurations and many more configuration interactions. 

At present X am limited to 61 even and 61 odd configurations and X try to 
include everything up through n * 9. The current program runs on Alpha 
workstations. X decided to test the big program on Fe I and Fe XX to see 
whether there was any great difference in the low configurations compared to 
those from the Cray program. 
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Besides increasing the number of El lines by a factor of 6 to 7.7 
million, there was an unexpected result; the electric quadrupole trams it ions 
were 10 times stronger than before because the transition integrals are 

weighted by r*2 they become very large for high n, and because there are 

numerous configuration interactions that mix the low and high configurations. 
As a check I was able to reproduce Garstang's (1962) lower results by running 
his three configurations with my program. Since my model atom is still only 
a subset of a real Fe II ion, the true quadrupole A values are probably 
larger than mine. The magnetic dipole lines are affected by the mixing but 
the overall scale does not change. 

Because of this discovery I decided that there was no point in computing 
the small array cases. I have been running with as many configurations as 
I can and with thousands of parameters in the Hamiltonian. The computer 
runs take much longer to set up and produce than I had expected. I have 
concentrated on redoing the low iron group spectra, especially to get data 
for supernova modelers. I have done only Ca I — Zn I, Ca II — Zn II, 

Cu I — Cu XXIX, Zn I - Zn XXX, for practice at high stages of ionization, 

C I, C II, S I, and Cl I and Ag I for people who were working on the 
laboratory spectra. Check my web site kurucz.harvard.edu for current 
additions. 

My latest calculations have been for carbon I and sulphur I, and 
silicon I is under way using the same elaborate approach as for Cl, 
which took many months to do. These line lists greatly increase the 
number of lines in the ultraviolet, in the visible, and especially in 
the infrared. They will increase the opacity in A, F, and G stars. 

They will account for many unidentified lines in the sun. Here are 
parts of the .README files: 

0600 .README 9 October 2004 0600 C I 

Carbon I has now been completely redone and expanded and put in 
directory /atoms/0600. The computation was divided into four sections: 
suffix X, roughly all levels up through n * 9-12 with all configuration 
interactions; suffix R for Ritz, all higher levels up though n « 20 for 
low 1 with no configuration interaction; suffix P for polarization, all 
higher levels up though n * 20 for high 1 with no configuration 
interaction; and, by extrapolation, levels up through n * 99 that 
connect to s2p2, sp3, 2p3s, 2p3p, 2p3d levels and merge into a continuum. 

All higher levels up to n « 20 were filled in by Ritz expansion or the 
polarization formula for average energies. The polarization average 
energies are also given for all configurations with 1<21 for all n<100. 

These are the configurations in each section. 
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20 October 2004 1600 S I 

The final results for the electric dipole, magnetic dipole, and 
electric quadrupole lines are in the files GFEMQ1600.POS for wavelengths 
in nm or in GFEMQ1600 .WNPOS for frequencies in wavenumbers in cm-1 for 
lines between observed levels. Only forbidden lines from metastable 
levels are included. There are 2161 even and 2270 odd levels with 
configurations up to n ■ 15 or 16. There are 225605 electric dipole 
lines. Of these 24722 lines are between known energy levels and have 
good wavelengths. The rest are predicted. 


From the hydrogen sequence data that have been published Z have been 
able to interpolate or extrapolate the energy levels and to generate the 
lines up to n = 99 for elements up through phosphorus. Isotopes are treated 
only through boron. The starting data for sulfur and higher are scant so 
that the interpolations and extrapolations may not be reliable. NIST should 
write a general program and compute all these data once and for all. 

I had an urgent request for hyperfine and isotopic splitting for Cu X so 
I incorporated the known splittings into my previously calculated Cu I line 
list and put the new data on the web. Fits to the solar spectrum look very 
good. 


I have copied my extant test spectrum calculations to the web 
directory. The spectra were calculated either at resolving power of 500000 
or 2000000 in Doppler space so that a shift in point number is a Doppler 
shift. Depending on effective temperature, the wavelengths range from 1 nm 
to 300$\mu$m. Some of the spectra have been broadened to lower resolution 
on the same point spacing and have been resampled to a smaller number of 
points. Eventually the resolutions will include 100000, 50000, 30000, 
20000, 10000, 5000, 3000, 2000, 1000, 500, 300, 200, and 100. Oftentimes 
users can find what they need "off the shelf" on the web site. As I fill 
requests for new calculations, I add the new spectra to the website. This 
is the current list: alpha Cen, Arc turns, A spectral types, Betelgeuse, 
Procyon, Sirius, Vega, GL411, Sun. 

Fiorella Castelli in Trieste and I used my new distribution function 
program to compute new opacity tables and then new grids of model 
atmospheres with improved convection. The grids of models, fluxes, and 
colors have been put on the Kurucz web site. 

We have had complaints that we should be modern and use the Opacity 
Project and Iron project opacities. We (Castelli and Kurucz 2004) did an 
experiment in which she installed the Iron Project Fe I photoionization 
and autoionization opacity from Bautista (1997) into my model atmosphere 
program ATLAS9. The opacity looks huge with many autoionizing resonances 
in the continuum that required a great deal of supercomputer time to 
produce. In the Kurucz line lists the autoionizing resonances are treated 
as line opacity on top of a smooth continuum opacity. I compute all the 
possible transitions between J levels while Bautista could treat only 
term-to-term transitions. In some intervals I have much more opacity 
because I have 39000 levels. For some terms Bautista has much higher 
(smooth) continuum opacity and mine are wrong. The physical temperature 
structures of models computed with the two opacity treatments are the same, 
which implies that the total opacity, where it matters, is the same 
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averaged over small wavelength intervals for Bautista and Kurucz. In 
general the flux is the same as well except at one small interval where 
Bautista is higher and another where Kurucz is higher. Agreement does not 
imply that the models are correct. My old line lists are incomplete (that 
is why 1 am recomputing them) , so there must also be opacity missing from 
Bautista, and there are still systematic errors in the energy distributions 
predicted from the models. Comparisons with the observations (Bell, 
Balachandran, and Bautista 2001; Cowley and Castelli 2002) indicate that 
both computations appear to understimate the real Fe I absorption. 

And by implication all the iron group spectra. 

I have done minimal travel during this grant because it would slow 
down the work on the line data. 
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